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ABSTRACT 

Wc present a near infrared study focused on the detection and characterization of the youngest 
stellar component of the NGC 604 giant star- forming region, in the Triangulum galaxy (M 33) . By 
means of color-color diagrams derived from the photometry of JHK S images taken with Gemini- 
NIRI, we have found 68 candidate massive young stellar objects. The spatial distribution of 
these sources matches the areas where previous studies suggested that star formation might be 
taking place, and the high spatial resolution of our deep NIRI imaging allows to pinpoint the 
star-forming knots. An analysis of the fraction of objects that show infrared excess suggests 
that the star formation is still active, supporting the presence of a second generation of stars 
being born, although the evidence for or against sequential star formation does not seem to be 
conclusive. 

Subject headings: M33 — Hn regions — individual: (NGC 604) — galaxies: starburst — stars: early-type 



1. INTRODUCTION 

In this paper we present a near infrared (NIR) 
photometric study of the NGC 604 star-forming 
region. This study is focused on the detec- 
tion and first characterization of the NGC 604 
newest generation of stars. NGC 604 is one of 
the most outstanding giant Hn regions (GHRs) in 
the Local Gr o up. W ith a Lhu = 2.6x10 39 ergs -1 
( Bosch et all l2002t ) it is the second most lumi- 
nous Hn region after 30 Doradus. The observation 
and study of these two regions play a fundamen- 
tal role in our understanding of the astrophysical 
processes associated with massive stars and, in 
particular, massive star formation and the envi- 



ronments where it takes place. At a distance of 
840 kpc, NGC 604 is located far enough to observe 
global characteristics but still close as to resolve 
individual objects and small structures (1"~4 
parsecs). NGC 604 and 30 Dor provide the link 
between Galactic GHRs, which we are able to 
study in full detail, and distant GHRs, whose 
global properties can be traced to cosmological 
distances. Although they share several common 
characteristics they also exhibit substantial dif- 
ferences such as their spatial structure and ages 
distribution. 

The ionizing stellar cluster is co mposed of a 
popu lation of at least 200 O-type stars (jHunter et al 
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19961 ). Wolf-Rayet (WR) and Of stars, or candi- 
dates, were first identified and classified in stud- 
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i es bylConti k Massevl (Il98ll) iD'Odorico k Rosa 



198lh . iRosa k D'Odoricol [19821), iDiaz et af 



19871 ). and iMassev k Jo 



hnson (1998). New WR 



objects were detected by iDrissen et al.l (|l993l ) by 
means of Hubble Space Telescope (HST) images, 
and an accu rate classificat i on of k nown WR stars 
was m ade in lDrissen et al.l ( 20081 ). iTerlevich et al.l 
(Il996l ) spectroscopically classified two objects, 
one of them is a transition star from Luminous 
Blue Variable (LBV) to WR and the other is 
a red supergiant (RSG) star. There are forty 
OB stars with accurate spectral classification by 
iBruhweiler et al.l (|2003h . three of them exhibit 
Of/WR signatures (wind profile UV emission lines 
of Nv, Civ, and Siiv togethe r with strong emis- 
sion line of He ii at A1640). lEldridge k Rclano 
(|201ll ) studied the known evolved stellar popula- 
tion in the region by means of Spectral Energy 
Distribution (SED) fitting of HST UVIJHK pho- 
tometry. The authors used the UVI photo met- 
ric measurements from iHunter et al.l (|1996l ) and 
they performed NIR photometry using archive im- 
ages from HST-NICMOS (NIC2) in the F110W, 
F160W, and F205W filters. The latter HST data 
set, was originally used by iBarba et al. ( 2009h 
to perform a photometric study of the region, 
in which they found five RSG candidates and 
a dozen candidate massive young stellar objects 
(MYSOs). However, the selection criteria for IR- 
excess sources was based on the maximum errors 
present in the NICMOS photometry. This sug- 
gests that a considerable number of sources may 
yet be detected with deeper images, using photo- 
metric errors appropriate for each object. 

The age of the main stell ar population was es - 
timated by several authors. Hunter et al. ( 19961 ) 
analyzed the color-magnitude diagrams obtained 
from UVI HST-WFPC2 photometry; the authors 
fitted isochrones and considered the age constrains 
imposed by the presence of WR candidate stars. 
They concluded that the average ages for the stars 
in NGC 604 range from 3 to 5 Myr, and that 
the presence of RSG could suggest the existence 
of an older subpopulation. This older subpop- 
ulation was later studied by lEldridge k Relanol 
(|2011l ) who concluded that the RSGs belong to 
a for mation episode that occurred 12.4 ±2.1 Myrs 



plorer (IUE), together with the nebular emission 
lines and the higher order terms of the Balmer 
series and He I lines in absorption (both obtained 
from optical spectra taken with the 4.2-m William 
Herschel Telescope). The authors applied three 
different techniques of evolutionary synthesis and 
photoionization models (optimized for young star- 
forming regions) and concluded, consistently by 
the three methods, that the main ionizing clus- 
ter has an age of ^3 Myr. The physical and 
ki nematical properties o f NGC 604 were studied 
in IMaiz- Apellaniz ( 2000h (and references therein), 
by means of spectroscopic data obtained with the 
WHT (in particular, the lines formed in the warm 
ionized gas) and images from HST-WFPC2. The 
author concluded that NGC 604 suffered one large 
starburst 3.0-3.5 Myr ago that formed the main 
ionizing cluster in the region. 

Regarding the stellar population spatial distri- 
bution in the region, it is worth mentioning that 
the objects in NGC 604 are not centrally con- 
centrated (as is the case for 30 Dor) but spread 
over a la rge projected area ( ~ 10000 pc 2 for clus- 
ter A in IHunter et al.l Il996l ) , composing a struc- 
ture called the Scaled OB A ssociation (SOBA) by 



Mai'z- Apellaniz et al.l (|2004f ) 



Acc ordin g to the studies performed bylTenorio-Tagle et al 
(|2000h and IMaiz- Apellaniz et ail (|2004l ). the in- 
terstellar medium exhibits two with differ- 
ent excitation state: a high excitation region sur- 
rounded by a low-excitation region, but not with 
a concentric geometry. These structures also show 
different kinematic behavior and physical condi- 
tions. On the west side (high ionization) there 
can be identified three cavities, which seem to be 
bubbles bursting into the M33 halo (in our di- 
rection). There is a boundary, often cited as 'the 
ridge', which extends along the north-south di- 
rection defining a clear east-west division within 
nebulae. The east side (low ionization) exhibits a 
quieter kinematic behavior. This dichotomy is also 



obser ved in the X-ray emission (jTullmann et al 



2008). 



Detaile d studies on t he dy namics were also per- 
formed by lYang et al.l (|l996l ). who found that the 
width of the integrated velocity profile needs the 
contributions of virial motion, thermal broaden- 



a S°- i Gonzalez Delgado k Perez| (|2000|) analyzed Wstellar winds and SNRs. I Medina Tanco e~l 



stellar wind resonance lines in ultraviolet spec 
tra taken with the International Ultraviolet Ex 



(|l997h performed a study on the dynamics of the 
region by means of the Ha emission and they 
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concluded that a few tens of WR stars would 
be enough to feed the necessary kinetic energy. 
They also stated that the turbulence is relatively 
young, giving another element in favor of recent 
processes of stellar formation taking place a few 
million years ago. Analyzing the velocity disper- 
sion versus intensity diag r am o ver the whole re- 
gion, unoz-Tunon et ail(|l996l ) found that most 
velocity dispersion values measured in NGC604 
are supersonic. The authors also used the same 
diagram as a tool to identify the zones correspond- 
ing to structures such as shells, loops or bubbles 
which are generated by the interaction of the mas- 
sive stars' winds with the interstellar medium. 

The N GC 604 molecular component was widely 
studie d bvlViallefond et al.l (Il992l). I Wilson fc Scoville 
19921) . lEngargiola et all (l2003l).lRosolowskv et al " 



pOOa ). and iTosaki et alj ^2007l)~ In the latter 



study the authors proposed a sequential star for- 
mation scenario triggered by the expansion of the 
Hil region. This idea was reinforced in a recent 
study with higher spatial resolution i n CO, HCN. 
and 8 9- GHz continuum emissions bv iMiura et al 
(|2010h . 



In this study we present new observations of 
NGC 604 acquired using the Gemini-North tele- 
scope with NIRI in the J, H, and K s hort {K s here- 
after) bands. The high quality data obtained al- 
lowed us to perform a photometric study to iden- 
tify and characterize candidate MYSOs in the ion- 
izing cluster , confi rming the candidates found by 
iBarba et al. I (120091) and increasing the sample pre- 
sented in that study by a factor of five. Our data 
and analysis methods arc described in Section [2] 
Results and implications are discussed in Section 

m 

2. DATA REDUCTION AND ANALY- 
SIS METHODS 

2.1. Observations 

The observations were made using the Gemini- 
North telescope with NIRI (Near Infrared Imaging 
and Spectrometer) during the nights of September 
08 and 11, 2005 (Proposal GN-2005B-Q-3). NIRI 
was used in the imaging mode with the f/6 cam- 
era. This configuration yields a field of view of 
120"xl20" with a resulting plate scale of 0.117 
arcsec pixel -1 . The images were taken under ex- 
cellent seeing conditions, averaging 0."35 FWHM 



measured in the J, H, and K s images. Tablc[T]lists 
the broad-band filters used, together with their 
central wavelengths and wavelength coverage in 
Columns 2 and 3, respectively; the exposure time 
for each individual image is given in Column 4, 
and in Columns 5, 6, and 7 we have included de- 
tails from the photometry that will be discussed 
in the following sections. 

The observing sequences were made following 
a square pattern with offsets of about 10" every 
two consecutive field images, interspersed with a 
sequence of sky images, also taken with dither- 
ing offsets. The dithering sequence allows the 
removal of bad pixels and other detector artifacts, 
although the effective area that gets fully exposed 
is smaller (107" x 107"). Figure [T] shows the color 
RGB image of NGC 604 composed from our final 
broad-band images: J (blue), H (green), and K s 
(red). 



2.2. Data Reduction 

Most of the NIRI images exhibit a vertical strip- 
ing pattern which originate in the array's elec- 
tronics. Although its spatial profile is well known 
(eight columns wide) , its intensity varies from one 
exposure to another, even within an individual im- 
age, and it differs between the four image quad- 
rants since each quadrant in the NIRI array is read 
independently. The routine nirinoise.py (kindly 
provided by Gemini Staff) was applied to each im- 
age (field images, sky images, flats on/off images, 
and darks) in order to remove this pattern. Al- 
though it was not possible to completely remove 
it in all individual images, the residual pattern, 
when it was present, was only of a few ADU and 
being a small spatial scale pattern, it was cor- 
rected in the final image, obtained by averaging 
the dithered images. The second correction was 
performed to account for the NIRI detector non- 
linearity. This correction must be applied to all 
images with count rates greater than ~ 20 ADU 
s . The routine, also supplied by Gemini Staff, 
was applied to all J, H, and K s field images. 
After these two initial image corrections the re- 
duction procedure that follows was pursued with 
IRAF0 routines in the gemini. niri package. Fol- 



X IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
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Tabic 1: Main characteristics of the broad-band filters and the photometry 



Filter 


A c (fim) 


Coverage (/im) 


Exp. (s) 


Mean mag_err 


Max. mag 


Min. mag 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


J 


1.25 


0.97-1.07 


120 


0.06 


23.0 


16.0 


H 


1.65 


1.49-1.78 


60 


0.07 


22.5 


15.5 


K s 


2.15 


1.99-2.30 


40 


0.10 


22.0 


15.5 



Note. — Column description: (1) broad-band filters used; (2) filter's central wavelength; (3) filter's wavelength coverage; (4) 
the exposure time for each individual image; (5) mean photometric uncertainties; (6) the photometric limiting magnitudes for 
an error smaller than 0.4 magnitude; (7) the magnitude at the saturation limit. 
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Fig. 1. — RGB composite image of NGC604 field: J (blue), H (green), and K s (red). North is up, East to 
the left. Field size at NGC604's distance is roughly 490x490 pc 2 . 
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lowing the instructions in the NIRI web pages, the 
images were sky-subtracted and flattened. Short 
dark exposures were used to generate a mask to 
identify bad pixels. Special care was taken when 
combining the individual images to obtain the fi- 
nal image with the IMCOADD task. Due to possible 
confusion generated by crowding and the highly 
variable background, the IMCOADD align method 
was set in 'user' and the script was called inter- 
actively twice: a first call to mark the reference 
objects in each image and generate the transfor- 
mation map, and a second call to generate all the 
transformed images as well as to derive the aver- 
age image, with bad pixels and flagged cosmic-ray 
hits omitted. The average image was the final im- 
age used to perform the photometry. 
Some images exhibited a spurious, large scale, 
light pattern, which in most cases increased the 
counts by no more than 0.5-1.0% above the av- 
eraged background counts in a extended area of 
the field. This effect seems to appear, every time, 
after the offset to observe the sky region. Being 
a low spatial frequency pattern it is not expected 
to affect our point-source photometry. Neverthe- 
less, we examined its effect by generating two final 
images for each filter, one including all the images 
and a second one including only those images with- 
out any trace of the illumination pattern (reducing 
the number of useful images by about 50%). We 
performed photometric analysis in both images 
and compared the results. As expected, the mea- 
sured magnitudes were the same for both images, 
confirming that they were not affected by the light 
pattern, but the stacked image made with fewer 
individual images turned out to be less deep and 
the measured magnitudes presented more scatter. 
We therefore decided to pursue the study using 
the photometry derived from the coadded image 
that included the images with the light pattern, 
discarding only a few in which the pattern was 
more evident. 

2.3. Photometry 

2.3.1. PSF Fitting 

As evident in Figure [Q the NGC 604 field is 
highly crowded and the background nebular emis- 
sion exhibits strong variations on small spatial 



scales. These conditions drove us to pursue Point 
Spread Function (PSF) photometry. 
The objects for PSF fitting were selected using 
the DAOFIND IRAF task from an image created 
by adding our J and K s images. These images 
were aligned using the GEOMAP and GEOTRAN 
IRAF tasks to derive and apply, respectively, the 
geometrical transformation (only linear shifts). 
By performing the object finding in the J+K s im- 
age we guaranteed the detection of stars in both 
extremes of the total wavelength range covered 
by the NIR filters. A threshold of five sky sigmas 
over the background value was set as the detection 
limit. 

Stellar magnitudes were obtained by PSF fitting 
over the selected stars using the ALLSTAR rou - 
tine from DAOPHOT software (IStetsonl . Il987l ) 



Universities for Research in Astronomy, Inc., under coop- 
erative agreement with the National Science Foundation. 



in IRAF. The analytic PSF model that gave the 
best fit was a Moffat25 function (FWHM-3.0 
pixels ^0."35) with parameters that could vary 
quadratically with position in the image. Al- 
though a standard procedure, PSF construction 
and fitting in crowded fields involves an iterative 
and careful process. As there are no isolated stars 
available to be chosen as bona fide PSF stars to 
construct the PSF model, it requires several se- 
quences of PSF modeling and neighbor star sub- 
traction with an analytic PSF model of higher 
order each time, before the final model can be ob- 
tained. Pixels used to determine the stellar profile 
were those with less than 12000 ADU as this is 
the highest value for which the linearity correc- 
tion is expected to be reliable. The few stars that 
exceed this limit in their central pixels had their 
magnitudes determined from fits using the wings 
of their profiles. The photometric uncertainties 
include the internal photometric error for each 
source as determined by DAOPHOT (consider- 
ing photon counts statistics, NIRI noise and gain 
characteristics and PSF- fitting errors), to which 
the aperture correction uncertainty and the error 
associated with the transformation to the Mauna 
Kea Observatory (MKO) standard system were 
quadratically added. Aperture correction was cal- 
culated for each filter by performing aperture pho- 
tometry (with the aperture used to measure the 
standard stars) over an average of ~ 30 PSF stars 
with their neighbors previously subtracted. The 
error in the aperture correction is ~ 0.02 mag. for 
the three filters. The uncertainties associated with 
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the transformation of instrumental magnitudes to 
the MKO standard system arc 0.08, 0.01, and 0.02 
mag. for J, H, and K s , respectively. In Table [T] we 
have listed, for each filter, the mean photometric 
uncertainties in Column 5, the photometric limit- 
ing magnitudes for an error smaller than 0.4 mag 
in Column 6, and in Column 7 the magnitude at 
the saturation limit. Figure [2] shows the error dis- 
tribution with magnitude for each filter. 



2.3.2. Transformation to the Standard System 

The instrumental magnitudes were transformed 
to the MKO infrared photometric system, by 
means of aperture photometry performed on the 
standard stars. As the observed standard stars 
were limited to only one standard star per filter 
used in each observing night, we looked for a way 
to check the zero points of our magnitudes. A di- 
rect comparison with our photometry and 2MASS 
photometry was not possible as the magnitude 
limits of 2MASS at SNR=10 are 15.8, 15.1, and 
14.3 for J, H, and K s filters, respectively; which 
correspond to the magnitudes of stars that sat- 
urate in our images (see Table [J). Nevertheless, 
a comparison was possibl e by m eans of a study 



published by ICioni et al. (|2008h . The authors 



performed wide-field JHK S observations of M33 
using WFCAM at UKIRT, referring the astro- 
metric and photometric calibration to the 2MASS 
point source catalog. Fortunately, there are 10 
sources in common with our photometry which 
allows us to compare their magnitudes with ours. 
The results of the comparisons are shown in Figure 
[3] for the J, H, and K s filters where the presence 
of a small magnitude offset between both sets can 
be seen. As no additional trend stood out as ev- 
ident, and since the source number is small, we 
have limited our comparison to an evaluation of 
zero-point offsets. Performing a simple average to 
estimate these offsets we found them to be 0.06, 
-0.01, and -0.13 mag for J, H, and K s , respectively. 
Although individually small, they add up to one 
tenth of a magnitude in the (H-K s ) color. These 
offsets were added to the zero-points calculated for 
our photometry. The errors on these zero-point 
offsets to the magnitudes of the MKO system are 
smaller than the photometric uncertainties. 




Fig. 2. — Error versus magnitude distribution for 
J, H, and K s photometry J band magnitudes 
show the smaller errors and K s the relatively larger 
ones. This Figure is presented in a color version in 
the electronic edition of the Astronomical Journal. 
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Fig. 3. — Comparison b etwee n our photometric 



results and ICioni et al.l (|2008|) for the set of 10 
stars in common. In the diagram there are plotted 
the magnitude differences for individual stars on 
each filter. Horizontal lines show the best fit for 
constant offsets at 0.06, -0.01, and -0.13 for J. H 
and K s , respectively. The error bar in the top 
left corner represents the mean uncertainty for the 
magnitude differences. This Figure is presented 
in a color version in the electronic edition of the 
Astronomical Journal. 
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2.3.3. Central and Field Regions. Photometric 
Uncertainties 

An analysis to examine the internal pho- 
tometric uncertainties calculated by ALLSTAR- 
DAOPHOT routine was made for each filter by 
means of adding artificial stars to the science im- 
ages. Our NGC 604 field was divided in two re- 
gions: the 'central region' and the 'field region'. 
As the main cluster of NGC 604 does not stand 
out as an evident increase in the stellar density, 
the central region limit was set using a NIRI im- 
age obtained with the Pa/3 narrow-band filter (de- 
tails on narrow-band images and analysis will be 
included in a forthcoming paper). This is illus- 
trated in Figure [4] where the smoothed Pa/3 con- 
tours at 5(7 used to define the central region limit 
are shown. The central region, enclosed by a ~150 
pc radius circle, is centered at a = 0l' l 34 m 33 s .14 
and S = +30°47'1."9 and its area (of 68000 pc 2 ) 
encompasses the NGC 604 SOBA. The region out- 
side the circle, the field region, has a surface of 
~118000 pc 2 and was used to account for field 
star contamination. 

A histogram of magnitude distribution (0.5 mag 
bin width) was generated separately for each re- 
gion. A new image was created by adding artificial 
stars to each region. The number of added objects 
represented 10% of the stars at each magnitude 
interval. The artificial star magnitudes were mea- 
sured following the same procedure employed for 
the natural stars. By comparing their measured 
magnitudes with their 'true' magnitudes, we found 
that the differences are in the range of the mag- 
nitude uncertainties calculated by the ALLSTAR 
routine. In Figure [5] we have plotted the differ- 
ence between the 'true' and measured magnitudes 
for the artificial stars for the J, H, and K s filters in 
the top, middle, and bottom panels, respectively. 
The bars represent the magnitude error calculated 
by ALLSTAR for the measured magnitudes. 

2.3.4- Limiting Magnitudes 

We have also used the magnitude histograms of 
the NGC 604 central region to define the faintest 
magnitude reachable for each band. The magni- 
tude limit was set equal to the magnitude of the 
first bin brighter than the peak of the histogram 
plot (e.g., see Figure 6). Following this criterion, 
the magnitude limits arc 21.5 mag for J band and 



20.5 mag for H and K s bands. Figure [6] shows the 
histogram obtained for the K s band where it can 
be observed that 20.5 mag bin comes right before 
the bin with the highest counts (21.0 mag). With 
these magnitude limits our sample is reduced to 
1649 objects. However, the photometric uncer- 
tainties of the remaining objects are reduced to 
0.04, 0.04, and 0.07 mag. for J, H, and K s , re- 
spectively (nearly half of the mean magnitude un- 
certainties for the whole sample, see Table [lj. 

2.4. Astrometry 

The astrometry was made in two steps, first 
an alignment in 'image coordinates' of the final 
J, H, and K s images; then a subsequent transfor- 
mation from image coordinate system to World 
Coordinate System (WCS). The transformation 
from 'image coordinates' to (a, S) celestial equa- 
torial coordinates was performed using CCMAP 
and CCSETWCS IRAF tasks, which computed the 
plate solution and created the WCS image, respec- 
tively. The astrometry was derived using 45 iso- 
lated stars, common to our images and version 
2.3.2 of the Guide Star Catalog (GSC). Special 
care was taken to select well separated stars in our 
NIRI field. The reference frame of the GSC cata- 
log is the International Celestial Ref e rence Frame 
and Equinox J2000.0 (|Lasker et all . l2008t l. The 
GCS 2.3.2 typical errors are 0."3 and the trans- 
formation RMS arc 0."35 and 0."25 for a and 6, 
respectively. 

3. RESULTS AND DISCUSSION 

The magnitudes obtained in the individual fil- 
ters were matched in a unique list containing 3627 
objects in the field in which all three J, H, and 
K s magnitudes were measured. Table [2] is an ex- 
cerpt of the photometry table available in its en- 
tirety as machine-readable table in the electronic 
edition of the Astronomical Journal. Column 1 
contains the internal object id; Columns 2 and 
3 the Celestial Equatorial coordinates (J2000.0) 
in degrees; Columns 4 and 5 the J magnitude 
and its uncertainty, respectively; Columns 6 and 
7 the H magnitude and its uncertainty, respec- 
tively; Columns 8 and 9 the K s magnitude and 
its uncertainty, respectively; and Columns 10 and 
11 list the (J-H) and (H-K s ) colors, respectively. 
The total effective area - taking into account the 
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dithering pattern - covered by our photometry is 
of 107" x 107" (-430x430 pc 2 at the distance of 
M33). 



3.1. Background/Foreground Stars 

As mentioned in Section 1, fortunately, NGC 604 
is situated in a privileged location from an obser- 
vational point of view: in a direction out of the 
Galactic plane and in the outer part of M33 (a 
galaxy which is orientated almost face-on); hence, 
it is expected that in the NGC 604 field there is 
little contamination with foreground and back- 
ground stars. We have estimated the contami- 
nation b y Galactic ob j ects i n our field using the 
model of iRobin et al. I (l2003lM which yields only 
six objects, most of them F dwarf stars belonging 
to the thick disk. Our main concern for the study 
presented in this paper is the contamination by 
field objects lie in the area of the color-color (CC) 
and color-magnitude (CM) diagrams occupied by 
MYSOs candidates, resulting in possible misiden- 
tifications and an overestimation of their number. 
To examine this effect we generated density plots 
(two dimensional histograms) for the CC diagram 
of the NGC 604 central region and the field re- 
gion surrounding it, as defined in the previous 
Section. Figure [7] exhibits the CC density maps 
for the central region (top panel), the field region 
scaled to match the central region area (middle 
panel), and the bin to bin subtraction resulting of 
the central region minus the field region (bottom 
panel). Each bin in these diagrams is a 0.05x0.05 
mag. square, in (J-H) and (H-K s ). There are 
three visible differences between the object distri- 
bution in the CC density maps from both regions, 
which become even more evident in the bottom 
panel: (a) the star population around (J-H)^ 0.0 
and (H-K s )^ 0.0, which is present in the diagram 
of the central region but missing in the field re- 
gion. These objects are located in the CC area 
belonging to early-type main-sequence (around 
06-B2 V) and early-type giants stars, hence they 
are expected to be part of the NGC 604 popula- 
tion and do not suffer noticeable reddening; (b) 
the number of objects that lie to the right of the 
reddening line for a 06-08 V star, they constitute 
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the IR-excess objects and are the main interest 
for this study. It is clear that the IR-excess ob- 
jects cover a larger area in the CC diagram of 
the NGC 604 central region than for the field re- 
gion, and the field shows almost no objects redder 
than (H-K s )~ 0.5; (c) the location of the red giant 
stars in the observed CC diagram, grouped around 
(J-fl)~1.0 and (H-K s )~0.25, changes due to the 
reddening introduced by NGC 604. This reflects 
in a relative over-density of reddened RGs when 
looking through NGC 604 and a relative absence 
of their unreddened counterparts, seen as negative 
counts in the circled region marked in the bottom 
panel. 



3.2. Objects with IR Excess: IR excess 
Fractions and Spatial Distribution 

IR excess in MYSOs arise from heated dust and 
gas located in the very young objects' surround- 
ings, being a circumstellar phenomena rather than 
interstellar. This material can be part of the origi- 
nal cloud material that was cither left behind after 
the protostellar collapse and not blown yet by the 
new born star or, it can be material which forms 
part of the accretion disk of th e forming object 
still in the contraction process (jStahler fc PallaL 
20041 ) . Therefore, the objects' location in the CC 
diagram allows us to detect protostellar object 
candidates, as sources that lie on the right side of 
the reddening line for a 06—08 V star, exhibit- 
ing an IR excess. Figures |8] and [9] show the CC 
and CM diagrams obtained with our photometry 
for the NGC 604 central region. Those diagrams 
include 693 selected objects within the upper mag- 
nitude limits at J=21.5, #=20.5, and K s =20.5 
mag. As there is an intrinsic (and variable) uncer- 
tainty in the determination of NIR colors for each 
object, we have considered as 'objects with IR ex- 
cess' those sources whose distance to the redden- 
ing line is larger than their uncertainties (at one 
sigma) in the (H-K s ) color, marked with filled cir- 
cles and squares in Figures [8] and [9] In particular, 
the square symbols in both plots denote objects 
that exhibit 'extreme IR excess', that is, their IR 
excess is larger than their errors by a factor of 
at least three. In the CM diagram we have also 
included, as a reference, the location of two well 
known Galactic objects: S1 06 IRS4, a Class pro- 
tostar (|Furuva et all 1 1999b and NGC 2024 IRS2, 
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Table 2: Final photometry table 



Id 


a(J2000.0) 


5(J2000.0) 


J 


o-j 


H 


o H 


K, 




(J-H) 


(H-K 8 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


2498 


23.6167607 


30.7957712 


21.79 


0.08 


21.18 


0.09 


21.42 


0.24 


0.61 


-0.24 


4171 


23.6168631 


30.7972395 


21.29 


0.08 


20.49 


0.06 


20.41 


0.12 


0.80 


0.08 


3628 


23.6169750 


30.7976702 


21.22 


0.07 


20.75 


0.08 


20.72 


0.16 


0.47 


0.03 


2387 


23.6170550 


30.7946332 


22.55 


0.14 


21.81 


0.17 


21.92 


0.38 


0.74 


-0.11 


2062 


23.6171246 


30.7913550 


21.58 


0.08 


21.15 


0.12 


21.03 


0.18 


0.43 


0.12 



Note. — Excerpt of the final photometry table available in its entirety as machine-readable table in the electronic edition of 
the Astronomical Journal. Column description: (1) internal object id; (2) and (3) Celestial Equatorial coordinates (J2000.0) in 
degrees; (4) J magnitude; (5) J magnitude uncertainty; (6) H magnitude; (7) H magnitude uncertainty; (8) K s magnitude; (9) 
K s magnitude uncertainty; (10) (J-H) color; and (11) (H-K a ) color. 




01 h 34 m 35 B 



30 s 



a (2000) 



Fig. 4. — Grayscale Pa/3 emission line image ob- 
tained with NIRI. The area of the central region 
in NGC 604 is indicated with a circle that encom- 
passes the contour at 5<r level Pa/3 emission, used 
to define the cluster radius. The area beyond the 
circle was used as the field region. 




16.5 17 17.5 18 18.5 19 19.5 20 20.5 21 21.5 
mag 

Fig. 5. — Results for J (top panel), H (mid- 
dle panel), and K s (bottom panel) bands, from 
ADDSTAR test to study the internal photometric 
errors. The bars represent the magnitude uncer- 
tainties calculated by ALLSTAR for the measured 
magnitudes of the artificial stars. This Figure is 
presented in a color version in the electronic edi- 
tion of the Astronomical Journal. 
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Central region 



Fig. 6. — Histogram plot of the K s magnitude dis- 
tribution in the NGC 604 central region. The lim- 
iting magnitude was derived from the histogram 
and it corresponds to the 20.5 mag bin, which 
comes right before the bin with the highest counts 
at 21.0 mag. 



originally identi fied as an Ultra- Compact Hll re- 
gion. Although iBik et all (|2003l) suggest that it 
might be a young massive star with a dusty disk, 
this type of object is still linked to very recent 
massive star formation. With the criterion pre- 
viously mentioned, we have selected a total of 88 
objects that exhibit IR excess, 32 of them showing 
extreme IR excess. In Table [3] there are listed the 
88 objects with IR excess, where Column 1 is the 
internal object id (as in Table [2]), Columns 2 and 
3 the Celestial Equatorial coordinates (J2000.0) in 
degrees, and Columns 10 and 11 list the ( J-H) and 
(H-Kg) colors, respectively. This table is available 
in its entirety as machine-readable table, in the 
electronic edition of the Astronomical Journal. 



Besides MYSOs, in GHRs we can also find 
massive evolved stars like WR stars, B[e] super- 
giants (sgB[c]) and Of stars. As was commented 
in Section 1, the evolved stellar population within 
NGC 604 has already been widely studied. Partic- 
ularly, the WR conte nt of NGC 604 was studied by 
Drissen et all (Il993h with the HST/WFPC. In a 
complete survey, the authors identified 14 Of/WR 
candidates (up to M B =22.88). From these 14 
Of/WR candidates we detected 11 objects, since 
in our NIRI images objects WR4a - WR4b and 
WR2a - WR2b could not be resolved, and WR11 




Field region 




Central region - Field region 
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Fig. 7. — Density maps in the color-color diagram 
for objects in the central region of NGC 604 (up- 
per panel) and the control field (middle panel). 
The bottom panel shows the bin to bin subtrac- 
tion resulting of the central region minus the field 
region. In this panel four features are evident, 
denoting the differences in the stellar populations 
of both regions: i) A relative over-density of early- 
type main-sequence objects around ( J-.ff)~0.0 and 
(H-K s )^0.0; ii) a relative over-density of red gi- 
ant stars (that belong to M33's disk population) 
with moderate extinction around (J-H)^1.0 and 
(H-K s )^0.25; iii) an area with negative counts, 
marked with a circle, due to the lack of unreddened 
red giant stars in the NGC 604 region; iv) the IR- 
excess objects that lie to the right of the redden- 
ing line for a 06-08 V star. These objects cover a 
larger area in the CC diagram of the NGC 604 cen- 
tral region than for the field region, and the field 
shows almost no objects redder than (H-K s )~0. 5. 



is under the detection limit (indeed it is also op- 
tically faint). From the 11 WR detected in our 
photometry, 7 show IR excess greater than three 
times their uncertainty in (H-K s ) color, whereas 
the remaining four do not exhibit IR excess at all. 
Therefore, from the 88 objects which exhibit IR 
excess, we have subtracted the 7 known WR stars. 

Regarding possible contamination with evolved 
B[e] stars, it is worth reviewing that the known 
objects with the B[e] phenomena are not abun- 
dant and most of the B [e] stars are relatively faint 
(jMiroshnichenkol 120061) . From the five classes of 
B [e] stars ( Lamers et all Il998l ) , we only expect to 
detect the most luminous objects (sgB[e]) at the 
distance of M33. Considering that the list of all 
known Galac tic sgB [e] candidates in the study by 
Krausl ( 2009h includes only 14 objects, we can infer 
that it is very unlikely to find a significant number 
of sgB [e] in NGC604. 

Summarizing, the NGC 604 central region con- 
tains 693 objects with 81 showing IR excess (with 
the known WR subtracted). There are 566 ob- 
jects in the field region (after the scaling cor- 
rection applied to match both surveyed areas), 
with 13 showing IR excess. Therefore, we can 
estimate a total number of 127 objects belong- 
ing to NGC 604 in our photometry. From the 
81 objects with IR excess in the region centered 
at NGC 604, we subtracted 13 objects to account 
the background/foreground contamination, which 
yields 68 objects with IR excess in the NGC 604 
central region, and a fraction of objects with JHK 
IR excess, rfm, of 68/127 (~ 54%). 

A s was previously m entioned in the Introduc- 
tion, Barba et al. ( 20091 ) found a dozen candidate 
MYSOs by means of HST-NICMOS NIR data. In 
the present study we were able to detect ten of 
these objects and increase the sample of MYSOs 
candidates by a factor of five. Before performing 
an object to object comparison it is important to 
bear in mind that the category of 'objects with IR 
excess' was assigned with different numerical cri- 
teria in the NIRI and NICMOS studies. While in 
the present study each object's IR excess was de- 
termined taking into account its uncertainty in the 
(H-K s ) color, in the NICMOS study a fixed limit 
of 0.4 mag was set as the minimum shift for all ob- 
jects to flag them as showing an IR excess. This 
limit was chosen considering the maximum uncer- 
tainties in the photometry of the objects included 



in NICMOS sample. Our comparison therefore 
focused on the 2 out of 12 objects selected as can- 
didate MYSOs in the NICMOS study that are not 
among the objects with IR excess in Table [3J The 
object located at a=23.64050, 5= +30.78196 lies 
extremely close to a RSG candidate which satu- 
rates in the NIRI images and might have ham- 
pered the chances of a clear detection in our pho- 
tometry. The object located at a=023. 63578, 5= 
+30.78483 is a faint source in a highly crowded 
area, two factors that affected the object detec- 
tion in NIRI images. An inspection of the colors 
for the 10 objects with IR excess in common in 
both samples shows that the objects exhibit dif- 
ferences particularly in the (J-H) color, mainly 
caused by a combination of dissimilar transmission 
curves between the NIRI- J and NICMOS-F110W 
filters, and the fact that several of these objects 
were close to the detection limit in the F110W fil- 
ter. 

Figure [10] shows the objects' distribution in 
the NGC 604 central region, where open circles 
are stars without IR excess and filled squares de- 
note the objects with IR excess. Contours corre- 
sp ond to the 8.44 GHz radio continuum emission 
bv IChurchwell fc Gosd (|l999t) . designated follow- 



ing the authors' notation and pink circles denote 
sub-regions for which rjin was calculated as it is de- 
scribed below. In FigurcHOlit is evident that most 
of the objects which present IR excess are located 
in the regions delineat ed by the radio con t inuum 
emission contours. In Churchwell &: Gosd (Il999h 
the authors concluded that the most luminous 
knots require the equivalent of five to eight 05 III 
stars to account for their ionization. An inspection 
of Figure [TU] yields a number ranging from 6 to 15 
embedded objects for each peak. The spatial dis- 
tribution of our candidate MY SOs is also in agree- 
ment with a previous study bv lTosaki etHI (|2007l ) 
who, by means of CO molecular emission observa- 
tions at CO (J = 1-0) and CO (.1 = 3-2), inferred 
that the region of the arc coincident with the radio 
continuum emission C-D and its extension to the 
East (see Figure 3 of the mentioned paper) have 
the appropriate high density and temperature con- 
ditions to favor and support ongoing massive star 
formation processes. In addition, for the region co- 
incident with the rad i o emi ssion peak at knot A, 
Maiz-Apellaniz et all (120041 ) . based on an exhaus- 



tive study of the different gas phases at NGC 604 
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Fig. 8. — Color-color diagram for stars in the central region of NGC604. Filled circles and squares denote 
objects with IR excess. Particularly, square symbols are objects that exhibit 'extreme IR excess'. As a 
visual aid, the intrins i c colo rs of unreddened stars according to their spectral type and luminosity classes 
taken from iTokun aga (2000) are also included. Straight lines indicate the reddening paths for the bluest 
and reddest stars, calcul ated considering an absorption of Ay = 10.0 mag and the extinction coefficients by 
Rieke fc Lebofskvl (|1985l ). Stars that lie to the right of the lower reddening line are those that show an 
infrared excess larger than one time and three times their photometric uncertainties, blue points and red 
squares, respectively (in the electronic color version of the paper). The error bars at ±1 sigma in the {J-H) 
and (H-K s ) colors were also plotted for each object. This Figure is presented in a color version in the 
electronic edition of the Astronomical Journal. 




Fig. 9. — Color-magnitude diagram for stars in the central region of NGC604. The same symbol codes as 
in Fig. [8] were used to flag stars that show IR excess. As a reference, we have added the expected location 
for stars accordi ng to their spectral types and luminosity classes, their (H-K s ) colors were taken the from 
iTokunagal (|2000l ). and we have used the magnitudes of Main-Sequ e nce st ars from lMartins fc P lez (2006|, the 
magnitude s of Zero Age M ain-Sequence stars from lHanson et aL (1997) and the magnitudes for supergiant 
stars from IWegnerl ( 2007 ). Straight line indicates t he reddening path c alcula ted considering an absorption 
of Ay=10.0 mag and the extinction coefficients by iRieke fc Lebofskvl (119851 ). We have also included two 
well known MYSOs, namely S106IRS24 and NGC2024IRS4 (see text for discussion). The error bars at ±1 
sigma in the (H-K s ) color and the K s magnitude were als o plotted for each obje ct. The distance modulus 



used (my-My = 24.62) mag corresponds to M33's distance (jFreedman et all 120011) . This Figure is presented 



in a color version in the electronic edition of the Astronomical Journal. 
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that these percentages were calculated with rel- 
atively few objects, and are therefore very sen- 
sitive to small changes in the magnitude limits 
or subregion sizes, the results just described are 
well established. Regions C, D, and B are located 
within the arc pointed-out bv lTosaki et al. ( 2007 ) 
to exhibit ong oing star-formation m e ntione d be- 



fore, whereas Maiz-Apellaniz et al. ( 2004 ) pro 



posed that region A may be witnessing induced 
star formation. On the other hand, region F, co - 
incidcnt with 'Cluster B' from Hunter et al. I (ll996l) 
has an age similar to the main central cluster in 
NGC 604, and E is coincident with the SNR found 



byE Odorico et al.l (|l98C ). If r\i r correlates with 
the age of the subregions, as is the ca se for studies 
in Galactic star-forming regions (see Haisch et al.l 
2000UHaisch et al.ll2001l ). the results of this analy- 
sis, based on the NIR photometric characteristics 
of individual objects, is consistent with the idea 
that most of the objects in regions A, B, C, and 
D belong to a younger generation than the main 
SOBA population, which is supported by the men- 
tioned previous studies based on global observa- 
tions of the gaseous component. 

3.3. Individual Objects 

We have also made a brief survey of the rele- 
vant individual objects present in our photometric 
study that were studied by other authors. These 
objects are identified in our NGC 604 field image 
in Figure [TT1 and plotted in the respective CC and 
CM diagrams in FigurefTp } Circles show WR stars 
included in the study of iDrissen et al.1 (|2008l ) and 
the references therein. Square boxes identify ob- 
jects resulted from a visual identification between 
our NIR I images and the UV i mages from the 
study by iBruhweiler et al.1 (|2003l ). These objects 
are mainly OB stars with spectral classification 
derived by the authors from UV HST-STIS spec- 
tra. The identification between the HST finding 
chart and our NIRI frames must be considered 
with the caveat that, at M33's distance, for some 
objects we might be looking at unresolved stellar 
blends. An illustrative example can be found when 
comparing with analog regions, s uch as the star 



formin g knots in 30 Doradus where IWalborn et al 
(1999) identified several strong knots lying 2-3 pc 
away from massive 03-6 stars. A similar con- 
figuration would fit within one arcsecond in our 
images. In Table [5] we have listed these objects, 



Table 3: Objects which exhibit IR excess. 



Id 


a(J2000.0) 


<5(J2000.0) 


(J-H) 


(H-K s ) 


(1) 


(2) 


(3) 


(4) 


(5) 


4234 


23.6323556 


30.7843609 


0.41 


0.72 


4063 


23.6325544 


30.7835649 


0.45 


0.46 


3302 


23.6326740 


30.7833725 


0.51 


0.96 


4361 


23.6328734 


30.7839901 


-0.05 


0.37 


4048 


23.6329417 


30.7830922 


0.47 


0.52 



Note. — 

Excerpt of the table listing objects which exhibit IR excess. 
The table is available in its entirety as machine-readable ta- 
ble in the electronic edition of the Astronomical Journal. 
Column description: (1) internal object id; (2) and (3) Ce- 
lestial Equatorial coordinates (J2000.0) in degrees units; 
(4) (J-H) color; and (5) (H-K s ) color. 



w o 
' — o 




a (2000) 

Fig. 10. — Spatial distribution of observed ob- 
jects in the NGC 604 central region. Radio con- 
tinuum contours and regions as identified by 



Churchwell fc Gossl (|l999h are included. Open cir- 
cles are stars without IR excess and filled squares 
denote the objects with IR excess. 
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where Column 1 is the Id. given in previous stud- 
ies preceded by a set of characters that identify 
the reference, Column 2 contains the object clas- 
sification from the previous studies, in Columns 
3, 4, and 5 we have included our measured J, -ff, 
and K s magnitudes, respectively. Those objects 
marked with c in Column 1 have a special com- 
ment in the following paragraph. 

Comments on individual objects: 

WR stars tend to be red objects in the NIR CC 
and CM diagrams since they are more luminous 
in the K band than the J band as they are being 
observed through their own obscured envelope and 
winds. 

D08-WR2 was first classified as WN bv lConti fc Massev 



Table 4: Infrared excesses derived for individual 
knots of star formation 

Region a(J2000) 5(J2000) tjir (excIR/all) 





C) 


C) 


(%) 


A 


23.639208 


30.782583 


69 (9/13) 


B 


23.639333 


30.780639 


75 (6/8) 


C 


23.634708 


30.782778 


75 (12/16) 


D 


23.633000 


30.783667 


88 (15/17) 


E 


23.637708 


30.777583 


(0/0) 


F 


23.632083 


30.788167 


(0/1) 


X 


23.635685 


30.785212 


15 (4/26) 



(|!981f ). who could not resolve it from D08-WR1. 
D08 did not include this object in their sample. 
B03 classified it as 04Iab-04Ia, since they could 
identify two objects at this position, being the 
most massive and luminous as well as the youngest 
objects of their sample. From a visual inspection 
of our images, it is evident that there is more than 
one object at this location although we are not 
able to resolve them photometrically. The (H-K s ) 
color of this object in our photometry is bluer 
than the expected for a WN object but also its 
K s magnitude is brighter than the magnitude of 
a single O I star, which can be explained as the 
composite emission from two (or more) unresolved 
objects, in agreement with B03 observations. 

D08-WR4 presents a similar situation as D08- 
WR2. This ob j ect w as classified as WN by 
Conti fc Massevi ( 198ll ). who could not resolve 
it from D08-WR3, and D08-WR5. D08 did not 
obtain a new spectrum for it, but at the same lo- 
cation B03 obtained a spectra and classified three 
objects: two are 0911 stars and a 091a. Again, 
in our images we are able to resolve only one ob- 
ject but a visual inspection suggests that at this 
location there is more than a single object. Be- 
sides, the blue shifted (H-K s ) color and the bright 
K s magnitude indicate that we are possibly mea- 
suring the composite emission of a tight group 
of massive stars (not necessarily evolved as WR 
objects) in agreement with B03 classification. 

D08-WR6 is one of the bright est objects in 
the re gion. It was first ob served by Drissen et all 
(|l993h and studied later bv lTerlevich et all (|l996f T 
who measured a He overabundance and extremely 



Note. — 



Each region is identified according to IChurchwell fe Gosj 
nomenclature, their center coordinates listed in 
Columns 2 and 3, their IR excess ratio shown in Column 4. 



in 

Q238 
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Fig. 11. — J band image with the identifica- 
tion of the stellar ob j ects i dentified in the liter- 
ature. Drissen et al. I (120081) sources are marked 
with square boxes and lBruhweiler et~ ( 2003 ) are 
shown with open circles. 
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Tabic 5: Stars with known spectral types 



Id 


Reference 


Classification 






J 








H 








K s 






3791 


D08-WR1 


WCE 


18 


.32 


± 





.03 


18 


.01 


± 





03 


17 


.00 


± 





,04 


1168 


D08-WR2 c 


WN 


17 


.00 


± 





.03 


17 


.01 


± 





.03 


1G 


.97 


± 





.04 


1168 


B03-867 A-B c 


04Iab-04Ia 
































3317 


D08-WR3 


WN 


17 


,68 


± 





,03 


17 


.52 


± 





,03 


17 


.19 


± 





.04 


3795 


D08-WR4 C 


WN 


16 


.49 


± 





,04 


10 


.35 


± 





,05 


1G 


.19 


± 





,06 


3795 


B03-578A-C C 


O9II-O9Ia-09II 
































1257 


D08-WR5 


WC6 


19 


.69 


± 





.03 


19 


.54 


± 





.05 


18 


.86 


± 





.06 


1265 


D08-WR6 


WNL 


17. 


.29 


± 





.02 


17 


.11 


± 





.02 


16 


.62 


± 





.03 


4021 


D08-WR7 


WC4 


19 


.26 


± 





.03 


19 


.12 


± 





.03 


18 


.78 


± 





.05 


4096 


D08-WR8 C 


WN6 


18 


.75 


± 





.03 


18 


.66 


± 





.03 


18 


.61 


± 





.05 


3255 


D08-WR10 


WN6 


18 


.85 


± 





.03 


18 


.79 


± 





.03 


18 


.68 


± 





.04 


4075 


D08-WR12 C 


WNE 


17 


.61 


± 





.03 


17 


.40 


± 





.03 


17 


.05 


± 





.04 


4075 


B03-690A-B C 


05111-BOIb 
































1247 


D08-WR13 C 


06.5Iaf 


18 


.40 


± 





.02 


18 


.45 


± 





.03 


18 


.48 


± 





.04 


1244 


D08-V1 


Of/WNL 


17 


,88 


± 





,02 


17 


.97 


± 





,03 


17 


.67 


± 





,04 


1631 


B03-117 


0411 


19 


,61 


± 





,02 


19 


.59 


± 





,03 


19 


.57 


± 





,06 


1545 


B03-238 


09III: 


20 


.67 


± 





,04 


20 


.45 


± 





,06 


20 


.30 


± 





.11 


1327 


B03-530A-E 


07V: 


18 


,37 


± 





,04 


18 


.43 


± 





,04 


18 


.42 


± 





,05 


4297 


B03-564 c 


0911 


18 


,34 


± 





.07 


18 


.28 


± 





.07 


18 


.11 


± 





.07 


4380 


B03-602B 


0911 


19 


.28 


± 





.11 


19 


.47 


± 





.08 


19 


.32 


± 





.10 


4362 


B03-675 


0711 


19 


.18 


± 





.04 


19 


.13 


± 





.06 


19 


.06 


± 





.10 


4224 


B03-899A 


0711 


18 


.73 


± 





.04 


18 


.02 


± 





.04 


18 


.40 


± 





.05 


4425 


B03-938 c 


06IIL: 


19 


.16 


± 





.09 


19 


.20 


± 





.11 


18 


.56 


± 





.13 


4335 


B03-951 c 


08V 


19 


.08 


± 





.10 


19 


.14 


± 





.11 


18 


.29 


± 
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Note. — Column 1 shows the object id. from this work, Column 2 includes different identifiers used in the literature for the 
same object, Column 3 lists the spectral classification and Columns 4 to 6 include N1R magnit udes derived in the present 
work. Those objects m arked with c in Column f, have a special comment in this Section. D08: Drisscn ct al. (2008), B03: 
Bruhwcilcr ct al. J2003T ), the authors have identified with same number and a different character those object with overlapping 
spectra. 
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Fig. 12. — Color-color (left panel) and color-magnitude (right panel) diagrams for stars with known spectral 
types. Reference tracks, reddening lines and known m assive young s t ellar objects are the same as those 
in Figures [8] and p Sq uares indicate objects listed by iDrissen et al.l ( 20081 ) and circles those studied by 

" Br 



Bruhweiler et al.l (|2003l) . This Figure is presented in a color version in the electronic edition of the Astro- 



nomical Journal. 



wide spectral emission lines (^2500 km s _1 ) as 
well as a spectral variability in a time span of <~ 10 
years. The authors classified this object as a tran- 
sition object between a Luminous Blue Variable 
and a WR star. 

D08-WR8 was already known as a WR object 
candidate when D08 precisely classified it as WN6. 
In the HST F170W images it is clear that there 
are two objects very close together in projection, 
although our photometry distinguishes only one 
object, that may cause the D08-WR8 blue color 
in our CC diagram. 

D08-WR12 is resolved into two objects in B03 
HST images, the authors classified them as 05III 
and a BOIb, although D08 is in disagreement with 
these spectral types, based mainly on the pres- 
ence of the Hell (A4686) emission line. D08 ar- 
gued that although their spectrum is a composite 
of two stars most of the emission lines must come 
from the brightest star (the BOIb for B03) and 
they classified it as a WN10 star. 

D08-WR13, this object is approximately 30 arc- 
seconds away from the central cluster. D08 sug- 
gested that its spectral type is 06 Iaf. This object 
does not exhibit IR excess in our photometry. 

B03-564 is really two stars but we are not able 
to resolve them, so our magnitudes have contribu- 
tions from both objects. 



B03-938 was classified as 06III:: by B03. In 
the HTS-STIS F170W image it is clear that there 
are two objects close together in projection while 
in our photometry only one is detected. Besides, 
these objects are located in a region within a 
ring of noticeable nebular emission which prob- 
ably causes the red color of this object in our CC 
diagram. 

B03-951, is located close to B03-938, in the 
same nebular ring, and it is also red for its 08 V 
spectral type according to B03 classification. 

4. SUMMARY AND CONCLUDING RE- 
MARKS 

We have performed a NIR photometric study 
on the NGC 604 massive star-forming region, using 
Gcmini-NIRI images. 

Our study was focused on the detection and 
individual identification of massive young stellar 
objects that comprise a new stellar generation in 
the region. 

By means of JHK S photometry we have iden- 
tified 68 candidates MYSOs. These objects were 
selected by their position in the CC diagram (each 
of these objects exhibits IR excess greater than 
their photometric error in (H-K s )). 

The analysis of the spatial distribution shows 
that most of the MYSOs candidates are located in 
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areas with strong nebular emission, delineated by 
the radio continuum contours at 8.44 GHz, where 
previous studies stated the possibility of star for- 
mation processes taking place. 

We have calculated the fraction of objects with 
IR excess throughout the whole region. The re- 
sult yields that several regions in NGC 604 show 
a high fraction of objects with IR excess for the 
cluster age. Although it is not possible to perform 
a direct comparison with star forming regions in 
our Galaxy, this could be interpreted as a recent 
burst of star formation. Studies of other GHRs are 
necessary to establish whether there is a general 
tendency for regions of massive star formation to 
present higher IR-excess fractions and to under- 
stand the true nature of IR-excesses in MYSOs. 

We have made a short review of individual ob- 
jects in NGC 604 included in our sample, which 
were already studied by other authors. 

All these results complete a general picture that 
exhibits the existence of a new stellar generation in 
the NGC 604 region, which is mainly taking place 
in the molecular gas that extends to the south-east 
side of the main SOBA population. 

This is the first study of the NGC 604 star form- 
ing region dedicated to the detection of individual 
MYSOs candidates. New detailed observations for 
each object are needed to confirm the nature of 
the MYSOs candidates we have found and also 
to provide information on the physical conditions 
of the circumstcllar matter which is causing the 
observed IR excess. 
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